Chapman-Enskog expansion and relativistic dissipative hydrodynamics 
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Starting from Boltzmann equation with relaxation time approximation for the collision term and 
using Chapman-Enskog expansion for distribution function close to equilibrium, we derive hydrody- 
namic evolution equations for the dissipative quantities directly from their definition. Although the 
form of the equations is identical to those obtained in traditional Israel-Stewart approaches employ- 
ing Grad's 14-moment approximation and second moment of Boltzmann equation, the coefficients 
obtained are different. In the case of one-dimensional scaling expansion, we demonstrate that our 
results are in better agreement with numerical solution of Boltzmann equation. 
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Relativistic fluid dynamics has been applied quite suc- 
cessfully to study and understand a wide range of col- 
lective phenomena observed in cosmology, astrophysics 
and the physics of high-energy, heavy-ion collisions. The 
earliest theoretical formulation of relativistic dissipative 
hydrodynamics also known as first-order theories (order 
of gradients), are due to Eckart [jj and Landau-Lifshitz 
Q. However these theories, collectively called relativis- 
tic Navier- Stokes (NS) theories, involve parabolic differ- 
ential equations and suffers from acausality and insta- 
bility. The Chapman-Enskog (CE) expansion has been 
the most common method to obtain first-order hydrody- 
namics from Boltzmann Equation (BE) Q. The deriva- 
tion of second-order fluid-dynamics by Israel and Stew- 
art (IS) from kinetic theory uses extended Grad's method 
[3| • The approach by Israel and Stewart may notguaran- 
tee stability but solves the acausality problem [5| at the 
cost of introducing two additional approximations: 14- 
moment approximation for the distribution function and 
use of second moment of BE to obtain evolution equa- 
tions for dissipative quantities. 

Grad's method, originally proposed for non-relativistic 
systems, was modified by Israel and Stewart so that it 
could be applicable in the relativistic case. In this ex- 
tension, known as 14-moment approximation, the distri- 
bution function is Taylor expanded in powers of four- 
momenta around its local equilibrium value. Truncating 
the Taylor expansion at second-order in momenta results 
in 14 unknowns that has to be determined to describe the 
distribution function. This expansion implicitly assumes 
a converging series in powers of momenta. In addition, 
it is assumed that the order of expansion in 14-moment 
approximation (expanded as a series in momenta) coin- 
cides with that of gradient expansion of hydrodynam- 
ics. This is evident because Grad's approximation trun- 
cated at second-order in momenta is not consistent with 
second-order hydrodynamics. 

Another assumption inherent in IS derivation is the 
choice of second moment of the BE to extract the equa- 
tion of motion for the dissipative quantities. This choice 
is arbitrary in the sense that once the distribution func- 



tion is specified, any moment of the BE will lead to a 
closed set of equations for the dissipative currents but 
with different transport coefficients. In fact, it has been 
pointed out in Rcf. [6] that instead of this ambiguous 
choice of the second-moment of BE by IS, the dissipative 
quantities can be obtained directly from their definition. 

In this communication, we present an alternate deriva- 
tion of hydrodynamic equations for dissipative quantities 
which does not make use of both these assumptions. We 
revisit the CE expansion of the distribution function us- 
ing BE in relaxation time approximation. Using this ex- 
pansion, we derive the first and second-order equations 
of motion for the dissipative quantities from their def- 
inition. For the case of one-dimensional scaling expan- 
sion, we compare our results with the results obtained by 
using both 14-moment approximation as well as second 
moment of BE (IS equations) . We demonstrate that our 
results are in better agreement with direct numerical so- 
lution of BE as compared to those obtained by using IS 
equations. 

Fluid dynamics is best described as a long-wavelength, 
low-frequency limit of an underlying microscopic theory. 
BE governs the evolution of single particle phase-space 
distribution function / = f(x,p) which provides a re- 
liably accurate description of the microscopic dynamics 
in the dilute limit. With this motivation, our starting 
point for the derivation of hydrodynamic equations is rel- 
ativistic BE with relaxation time approximation for the 
collision term Q 



(1) 



where, p^ is the particle four-momentum, is the fluid 
four-velocity and tr is the relaxation time. We define 
the scalar product u ■ p = u^p^. With / —> f and 
fo —> fo, Eq. (JTJ) describes the evolution of distribution 
function for antiparticles /. The equilibrium distribu- 
tion functions for Fermi, Bose, and Boltzmann particles 
(r = 1,-1, 0), is fo = [exp(/3 u ■ p — a) + r]^ 1 and for an- 
tiparticles is fo — [exp(/3 u ■ p + a) + r]^ 1 . Here, (3 = 1/T 
is the inverse temperature and a = fj,/T is the ratio of 
chemical potential to temperature. 
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In the CE expansion, the particle distribution func- 
tion is expanded about its equilibrium value in powers of 
space-time gradients. 

f = fo + Sf, Sf = <5/« + 6fW + . . . , ( 2 ) 

where Sf^ is first-order in gradients, 5f^ is second- 
order, etc. The Boltzmann equation, ([T]), in the form 
/ = /o — (Tj^/u-p)p^d f _ l f, can be solved iteratively as [1] 



fi — fo~ 



u ■ p 

where /i = /o + Sf 1 - 1 



P^d^fo, h = fo- 



tr 



u ■ p 

and h = fo + Sf w + 5f( 2 
first and second-order in gradients, we obtain 

Sf (1) = - — p»dJo, 



Sf 



(2) 



U ■ p 
tr 



p»p v 



dj — —d v fr 
\u ■ v 



(3) 
To 

(4) 
(5) 



u ■ p \u- p 

The above treatment to obtain Sf is valid for 8f (an- 
tiparticles) as well. 

For the sake of comparison, we also write down the 
Grad's 14-moment expansion [9( in orders of momenta as 
suggested by IS Q in orthogonal basis , 

Sf = /o/o (A n n + X n n a p a + \^ al 3P a p P ) + 0(p 3 ), (6) 

where, fo = 1 — r/o and An, A n , are determined from 
the definition of the dissipative quantities, Eqs. (p?|)-(ITTl). 
Since hydrodynamics involves expansion in orders of gra- 
dients, hence for consistency, CE should be preferred over 
14-moment approximation in derivation of hydrodynamic 
equations. 

The conserved energy-momentum tensor and particle 
current can be expressed in terms of distribution function 

as 



11] 



T< Mt = dpp»p v (f + f) = 



dpp"(f ~ f) =nu» + n", 



(7) 



where dp — gdp/[(2n) 3 y/^ + rn?], g and m being the 
degeneracy factor and particle rest mass. In the tensor 
decompositions, e, P, n are respectively energy density, 
pressure, net number density, and = g^ v — u^u v 

is the projection operator on the three-space orthogo- 
nal to the hydrodynamic four- velocity u^ defined in the 
Landau frame: T^ v u v — ew M . The metric tensor is 



fj 



flV — 



diag(+, — , — , — ). The bulk viscous pressure (II), 
shear stress tensor (w^") and particle diffusion current 
(n p ) are the dissipative quantities. 

Energy-momentum conservation, d fl T fJ ' 1 ' = and cur- 
rent conservation, d^N^ = 0, yields the fundamental evo- 
lution equations for n, e and u M 



e+(e + P + II)0- tt^V 



0, 



(e + P + Tl)u a - V a (P 



n) 



A^TT^ = 0, 

i9 + dpTf = 0. 



(8) 



We use the standard notation A — u^d^A for co- 
moving derivative, V" = A^ a 9 M for space-like deriva- 
tive, 6 — d^u^ for expansion scalar and A^ a B^ — 
{A a B f} + A /} B a )/2 for symmetrization. 

Even if the equation of state relating e and P is pro- 
vided, the system of Eqs. ([5]) is not closed unless the dis- 
sipative quantities II, n M and are specified. To obtain 
the expressions for these dissipative quantities, we write 
them using Eq. in terms of away from equilibrium 
part of the distribution functions (Sf, Sf) as 



n 



dpp a pP(Sf + Sf) , 



dpp a {5f-Sf) , 



tt^ = A^ / dpp a pP(5f + Sf) 



(9) 
(10) 



(11) 



where A^ = [A£A£ + A^A' - (2/3)A^A aj9 ]/2. 

The first-order dissipative equations can be obtained 
from Eqs. (|9>(fTTj) using Sf = from Eq. g]) 



n 



A 



a/3 



K dpp a 



i!pp a p 13 

TR 



P^ (f + fo) 



AU.U 



u.p 

dp p a pP 



u.p 

P^ (f - fo) 
tr 



u.p 



P^d^ (f + fo) 



(12) 
(13) 
(14) 



Assuming the relaxation time tr to be independent of 
four-momenta, the integrals in Eqs. (fl"2"1) - (|14l) reduce to 



n = -T R f3 n 9, n" = TRpnVa, -nT = Itr^g^ , (15) 

where a^ u = A^V Q it' 3 . The coefficients /3n, Pn and (3^ 
are found to be 



/Sn 

A, 

/9W 



i(l-3c^)(e + P)-H( e -3P) 



P(e + P) 
4P e - 3P 
~5~ + 15 



2(1)„ m* 2 
3/3 + 3/3 {[ P) '°- 



15 



((u-p)- 2 ) 



o+ 



(16) 
(17) 

(18) 



where (• • • ) 0± = / dp(- ■ ■ )(/ ± / ), and c 2 s = (dP/de) s/n 
is the velocity of sound squared (s being the entropy den- 
sity). It is interesting to note that these coefficients are in 
perfect agreement with those obtained in the Ref. 0] in 
which the evolution equations arc derived directly from 
their definition. This is due to the fact that in Ref. [f|, 
the coefficients /3n, /3„ and /3 ff , are associated with first- 
order terms and do not involve 14-moment approxima- 
tion. In the massless limit, ^ = rj/r^ = 4P/5 is also in 
agreement with that obtained in Ref. @ employing CE 
expansion in BE with medium-dependent masses. 
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To obtain second-order equations, we find that CE ex- 
pansion for the distribution function does not support 
derivation of hydrodynamic evolution equations from ar- 
bitrary moment choice of BE. Using the definition of 
dissipative quantities to obtain their evolution equa- 
tions comes naturally when employing CE expansion as 
demonstrated while deriving first-order equations (|15p. 
Second-order evolution equations can be obtained di- 
rectly by substituting Sf — + Sf^ from ((4| and 
([5|) in Eqs. ([9"T)-(fTTj). However, we follow the methodol- 
ogy discussed in Ref. |6j to derive second-order equations 
for the sake of convenience and clarity which will become 
evident in the following. We have confirmed that both 
these methods lead to the same result. 

We express the evolution of the dissipative quantities 
using Eqs. (|5])-(HT]) as 



n 



A a/ 3 



d PP a pP(Sf + 5f), 



<"> = K / dp P a (Sf-sf), 



(y.v) 



dp P a pP(Sf + Sf), 



(19) 
(20) 



(21) 



where the notations A^> = A^A V and B<^> = A^B"? 
represents space-like and traceless symmetric projections 
respectively orthogonal to u^. 

The co-moving derivative of the away from equilibrium 
part of the distribution function (Sf) can be obtained by 
rewriting the BE (JlJ in the form 

5f = -fo-— P^-yf-—, (22) 

U.p Tr 

and similarly for antiparticlcs. Using this expression for 
Sf in the second-order evolution equations ([T9l) - (f2"T]) . we 

get 



n 



n 

TR 



A 



a/3 



3 

-(fo 



dpp a p^ 



(fo + - 

V U.I 



fo, f /) 



(23) 



TV 
TR 



K d PP a (fo + — P 7 V 7 / 
J LV u.p 



(fo ->fo, f~> f) 



(24) 



TR 



fo 



+ (fo -> fa, f -> /) 



u.p 



-P 7 V 7 / 



(25) 



It is clear that in Eqs. (j23]) - ([25j) . the Boltzmann relax- 
ation time tr can be replaced by the relaxation times 
for the individual dissipative quantities Tn, T n , T n . At 
this stage, it seems as though the three relaxation times 
Tn, T n , TV are all equal to tr. This is because the col- 
lision term in the BE (JT|) is written in relaxation time 



approximation which does not entirely captures the mi- 
croscopic interactions. This apparent equality vanishes if 
the first-order equation (H"5j) is compared with the rela- 
tivistic Navier-Stokes equations for dissipative quantities 
(n = -(<?, = XTWa and tt^ = 2rja^). The dis- 
sipative relaxation times are then obtained in terms of 
first-order transport coefficients £, A, rj which can be 
calculated independently taking in account the full mi- 
croscopic behavior of the system 1^, HI ■ 

For the dissipative equations to be second-order in 
gradients the distribution function in right hand side of 
Eqs. (| 2 3 1) - ([25]) need to be computed only till first-order 
(Sf = Sf^) via CE expansion. We have repeatedly 
used the first-order equations for the dissipative quan- 
tities, Eq. (fT5|). to replace -> II, Va — > and 
a^ v — > 7r' Il/ such that the relaxation times appearing on 
the right hand side of Eqs. (|23l - ([25)) are absorbed. 

The second-order evolution equations of the dissipative 
quantities are finally obtained as 



n 



,w 



Put 



n 

T Un n ■ u 
— + /3„V"a 

Tn 

A„ n nv"a - ; 



<5nnff6> + ^UnT:^ 



\u n n ■ Va - 



n v uj 



£n n d ■ n , 



"V,.a 



' 717T n v 



A^B it 1 



(26) 



(27) 



r (v v ) — 



+ 2/3^(7 

Tir 



fMU 



'7T7T n y U 



- S^n^e + X^nUa^ - t. 



(28) 



We remark that although the form of the evolution equa- 
tions for dissipative quantities obtained here, Eqs. I[2"6"|)- 
(|2"5|) , are the same as those obtained in the previous calcu- 
lations using both 14-moment approximation and second 
moment of BE 14j . the coefficients obtained are different. 
In the following, we refer to the results in Ref. 14| as 



the IS results although the power counting scheme differs 
from the one employed originally by IS. 

For the special case of a system consisting of single 
species of massless Boltzmann gas, we find that 



=~ T 
Hit - , 1 ix 

5 



10 4 
7' ^- 3 ; 



(29) 



whereas these coefficients as obtained via IS approach are 



B IS = - 



IP 



= 2, 



<5 IS = 



(30) 



In this limit, although the coefficient of ir^O are same 
for both the cases (<5 W7r = d^ v ), the coefficient of cr A " / and 

tt#V>t are different ^ $f , ± t®). 

We note that CE expansion, as opposed to 14-moment 
approximation, can be done iteratively to arbitrarily 
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FIG. 1: (Color online) Time evolution of P L /P T in BAMPS 
(dots), IS (dashed lines), and the present work (solid lines), 
for isotropic initial pressure configuration (no = 0). 



higher orders. Hence using CE expansion, dissipative 
hydrodynamic equations of any order can in principle be 
derived. To obtain nth-order evolution equations for dis- 
sipative quantities, Sf = 5/W + 5f & + ■■■ + tf/^" 1 ) 
should be used in Eqs. ([2"5 |) -([2"5 ]l . For instance, substi- 
tution of Sf = 5/W + 5fW from Eqs. @ and © in 
Eqs. p3]l -([25 ]l will eventually lead to third-order evolu- 
tion equations. Derivation of third-order hydrodynamics 
as outlined above is left for future work. 

To demonstrate the numerical significance of the new 
coefficients derived here, we consider evolution in the 
boost invariant Bjorken case of a massless Boltzmann gas 
(e = 3P) at vanishing net baryon number density [l5j |. In 
terms of the co-ordinates (r, x, y, 77), where r = \/t 2 — z 1 
and r\ = tanhT 1 (z/t), the initial four- velocity becomes 
u' 1 = (1,0,0,0). For this scenario, II = = 0, and the 
evolution equations for e, 7r = —T 2 ir rin reduces to 

de 1 



(It 
dn 

~dr 



T 

n 



(e + P-n) 



3r 



The second-order transport coefficients simplifies to 



1 _ 38 



A IS = 2. 



(31) 
(32) 

(33) 



Initial temperature Tq = 500 Mev at proper time 
To = 0.4 fm/c are chosen to solve the coupled differential 
Eqs. (|3"Tj) and (|32l) . These values corresponds to LHC 
initial conditions [16|. We assume isotropic initial pres- 
sure configuration i.e. ttq — 0. Fig. [TJ shows the proper 
time dependence of pressure anisotropy which is defined 
as 

P L (r) P(t)-tt(t) 



The solid and dashed lines represents our results and the 
results from IS theory, respectively. The dots corresponds 
to the results of a transport model, the parton cascade 
BAMPS [HI El. The calculations in BAMPS are per- 
formed with variable values for the cross section such that 
the shear viscosity to entropy density ratio is a constant. 

We note that the results from IS theory always over- 
estimates the pressure anisotropy as compared to the 
BAMPS results even for viscosities as small as ij/s = 
0.05. The figure also shows that our results are in bet- 
ter agreement to the transport results. For very high 
viscosity, r)/s = 3.0, although at early times we have a 
better agreement with BAMPS as compared to IS, at 
later times there is a large deviation. This disagreement 
may be due to the relaxation time approximation for the 
collision term considered in Eq. (JTJ. 

Relaxation time approximation for the collision term 
assumes that the effect of the collisions is to restore the 
distribution function to its local equilibrium value expo- 
nentially. This is a very good approximation as long as 
the deviations from local equilibrium is small. For large 
•q/s values, shear viscous pressure is also large. This in 
turn implies that the deviations from equilibrium is not 
small resulting in breakdown of the relaxation time ap- 
proximation. It is important to note that large values of 
■q/s (> 0.4) are not relevant to the physics of strongly 
coupled systems like Quark Gluon Plasma (QGP). The 
QGP formed at RHIC and LHC behaves as a near perfect 
fluid with a small estimated rj/s « 0.08-0.2 [Ulifl]. The 
fact that we get reasonably good agreement with BAMPS 
results for rj/s < 0.4 (Fig. Q} also suggests that BE with 
relaxation time approximation for the collision term can 
be applied in understanding the hydrodynamic behaviour 
of QGP formed in relativistic heavy-ion collisions. 

To summarize, we have presented a new derivation of 
relativistic second-order hydrodynamics from BE. We use 
Chapman-Enskog expansion for out of equilibrium dis- 
tribution function instead of 14-moment approximation 
and derive evolution equations for dissipative quantities 
directly from their definitions rather than employing sec- 
ond moment of Boltzmann equation. In this new ap- 
proach, we get rid of two powerful assumptions of IS 
kind of derivation which is 14-moment approximation 
and choice of second moment of BE. Although the form 
of the evolution equations remain the same, the coeffi- 
cients are found to be different. Our results show a bet- 
ter agreement with the parton cascade BAMPS for the 
Pl/Pt evolution. 
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